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ABSTRACT 


Radiant  heat  transfer  was  analyzed  for  tents  consisting  of  single 
layer  fabric  walls  which  are  capable  of  partially  transmitting  thermal 
radiation.  External  radiant  temperatures  were  uniform.  Radiant  heat  loss 
from  the  tent  surface  to  the  external  surroundings  was  found  theoretically 
to  represent  approximately  25$  of  the  total  heat  loss  from  the  tent. 
Theoretical  predictions  of  the  fraction  radiated  to  the  floor  of  the  tent, 
based  on  experimentally  determined  surface  temperatures ,  agreed  to  within 
an  order  of  magnitude  with  the  measured  radiant  heat  flux  to  the  floor. 
Signal  noise  and  measurement  precision  affected  these  experimental  results. 
An  empirical  equation  for  predicting  sky  temperature  was  compared  with 
measured  sky  temperature  and  was  found  to  underestimate  the  sky  temperature 
by  approximately  10$. 


RESUME 


On  a  analyst  le  transfert  de  chaleur  rayonnante  de  tentes 
comportant  une  seule  toile  capable  de  transmettre  partiellement  le 
rayonnement  thermique.  Les  temperatures  de  rayonnement  ext^rieures  etaient 
uniformes.  La  provision  th£orique  du  transfert  de  chaleur  rayonnante  k 
travers  le  tapis  de  sol  de  la  tente,  bas4e  sur  des  temperatures 
superf icielles  mesurees  experimentalement ,  concorde  avec  le  flux  de  chaleur 
rayonnante  mesure  sur  le  tapis  de  sol.  A  cause  du  bruit  present  dans  le 
signal  et  de  la  precision  insuffisante  des  mesures ,  on  a  pu  verifier 
seulement  que  les  resultats  de  l'analyse  theorique  et  les  mesures 
exper imentales  etaient  du  meme  ordre  de  grandeur.  On  a  etabli  que  la  perte 
de  chaleur  rayonnante  vers  l'exterieur,  k  la  surface  de  la  tente, 
representait  environ  25$  de  la  perte  de  chaleur  totale.  Les  resultats 
obtenus  &  l'aide  d’une  equation  empirique  permettant  de  calculer  la 
temperature  du  rayonnement  ambiant  ont  ete  compares  avec  des  mesures  de 
cette  temperature.  Liquation  empirique  sous-estime  la  temperature  du 
rayonnement  ambiant  d'environ  10$. 
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GLOSSARY 


a,b  disk  radii 

A  surface  area,  m2 

F-.  view  factor  from  surface  i  to  surface  j 

h1J  vertical  separation  between  concentric  disks 

r..  a  vector  between  points  on  surface  i  and  surface  j 

R1J  radiative  heat  flux,  W/m2 

S  a  contour  line  around  a  surface 

T  absolute  temperatue,  K 

qg  radiative  heat  flux  of  the  pyranometer  sensor,  W/m2 

qnet  total  net  radiative  heat  flux,  W/m2 

Sji  Kronecker  delta 

e  J  surface  emissivity 

<j>  angular  step  size  around  the  vertical  axis  disk,  radians 

p  surface  reflectivity 

o  Stefan-Boltzmann  constant  =  5.67x1 0“ e  W/m2K1* 

x  surface  transmissivity 

0  angle  between  the  surface  normal  and  the  vector  r^j 


Subscripts 


A  ambient 

dp  dew  point 

i,j  surface  identifiers 

I  inside 

left 
outside 
right 
s  sensor 
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1  .0  INTRODUCTION 


This  study  examines  the  exchange  of  radiant  thermal  energy 
between  surfaces  of  a  tent.  A  general  theoretical  analysis  of  radiant  heat 
transfer  is  outlined  with  descriptions  of  techniques  required  to  perform 
the  required  analysis.  The  difference  between  this  problem  and  most  other 
radiant  heat  transfer  problems  is  that  the  fabric  surfaces  of  the  tent 
transmit  some  portion  of  the  incident  thermal  radiation  to  or  from  the 
external  surroundings  whereas  typical  engineering  materials  have  no 
transmission  of  radiant  thermal  energy.  Therefore,  the  conventional 
radiant  heat  transfer  analysis  [1]  has  been  extended  to  include  the 
transmission  of  radiative  thermal  energy. 

To  illustrate  the  technique,  a  conical  tent  shape  was  assumed. 

The  dimensions  of  the  theoretical  tent  were  chosen  such  that  the 
theoretical  tent  approximated  a  Canadian  Forces  5-Man  Arctic  Tent  as 
closely  as  possible.  For  experimental  verification  of  the  analysis,  a 
5-Man  Arctic  Tent  was  instrumented  to  measure  surface  temperatures  and 
radiant  heat  transfer  at  the  floor. 

In  this  study,  it  is  assumed  that  the  external  radiant 
temperature  is  uniform.  This  precludes  solar  heating  as  well  as  cases 
where  sky  temperatures  vary  substantially  from  the  ground  temperature. 

These  cases  are  to  be  studied  in  a  subsequent  investigation. 

Appendix  A  provides  interested  readers  with  the  computer  program 
used  to  do  the  numerical  computations  involved  in  determining  the  radiative 
heat  fluxes  in  the  test.  Appendix  B  gives  some  examples  of  the  viewfactor 
algebra  required  for  implementing  the  computer  analysis.  Appendix  C  gives 
a  listing  of  several  standard  library  routines  which  are  called  by  the 
program  in  Appendix  A  for  interested  readers  who  may  not  have  them  on  their 
computer  system. 

This  study  is  part  of  a  larger  investigation,  the  goals  of  which 
are  to  produce  an  understanding  of  the  important  heat  transfer  mechanisms 
in  tentage  and  to  provide  a  means  of  analysing  the  heat  transfer  in  tents. 
This  information  could  be  used  to  aid  in  the  development  of  tents  with 
superior  performance  characteristics  which  are  required  to  meet  the  ever- 
expanding  requirements  of  the  military  and  civilian  markets. 


2.0  THEORY 


The  equations  governing  radiant  heat  transfer  between  surfaces 
are  well  established  for  opaque  surfaces  [l,2]  but  seldom  is  transmission 
included  in  the  analysis.  For  most  engineering  materials,  the  transmission 
is  zero,  however,  surfaces  made  of  textile  fabrics  can  have  trans¬ 
missivities  between  0.05  and  0.25  [ 3 ] •  The  analysis  outlined  below  follows 
conventional  approaches  but  has  been  extended  to  include  the  transmission 
component  of  thermal  radiation. 


2.1  RADIANT  HEAT  TRANSFER  EQUATIONS 


The  analysis  begins  by  assuming  an  arbitrary  enclosure  as  shown 
in  Figure  1.  Let  RLIi  be  the  radiative  heat  flux  leaving  surface  i  and 
entering  the  enclosure,  RRIi  be  the  radiant  heat  flux  incident  on  the  inner 
surface  of  surface  i,  RRQi  be  the  radiant  heat  flux  leaving  surface  i  to 
the  external  surroundings  of  the  tent,  and  R^oi  l3e  the  radiant  heat  flux 
arriving  at  exterior  of  surface  i  from  the  surroundings. 

The  radiant  energy  arriving  at  the  exterior  side  of  surface  1  is 

given  by: 


(1 ) 


where  it  is  assumed  that  the  emissivity  of  the  external  surroundings  of  the 
enclosure  is  1.0  and  that  T^  is  the  radiant  temperature  of  the  external 
surroundings.  The  radiant  energy  leaving  the  exterior  side  of  surface  i  is 
due  to  emitted,  reflected  and  transmitted  radiant  energy  and  is  given* by: 
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Figure  1 .  Configuration  for  the  development  of  the  radiant  heat 
transfer  equations  for  an  enclosure  with  transmitting 

walls . 
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Similarly,  the  radiant  energy  leaving  the  interior  side  of  surface  i  is 
given  by: 


RLIi 


LoTi  +  piRRIi  +  t’’R 


i  LO 


(3) 


The  radiant  energy  arriving  at  the  interior  edge  of  surface  i  is  the  sum  of 
all  the  radiant  fluxes  leaving  the  j  internal  surfaces  of  the  enclosure 
which  strike  surface  i  directly.  Using  view  factor  algebra,  it  can  be 
shown  [1]  that  this  may  be  expressed  as: 


RRIi 


?  Fij  RLIj 


(4) 


Using  equations  1,3  and  4,  all  but  one  unknown  radiative  flux,  R^Ij '  can  be 
eliminated  from  a  system  of  simultaneous  equations  leaving: 


Z  (6tj  -  PiFi j )  RLij  =  eicrTi  +  TioTA 


(5) 


where  6^  is  the  Kronecker  delta  which  assumes  the  value  1  when  i  equals  j 
and  zero  otherwise . 

*  Solving  equation  5  for  RLj,  allows  the  remaining  radiant  fluxes 

to  be  evaluated  from  the  above  equations.  The  net  heat  loss  to  the 
exterior  surroundings  and  the  net  heat  loss  from  each  surface  are  given 
by: 


ROUTi”  RROi  “  RLOi 


(6) 


QNETi  ’  (RROi  “  RL0i)  +  (RLIi  “  RRIi) 


(7) 


Equation  5  can  be  solved  by  matrix  inversion,  Gaussian 
elimination!  or  by  iterative  routines.  Gaussian  elimination  is  the  mqst 
common  technique  for  small  systems  of  linear-  equations  and  was  the 
technique  used  here.  Appendix  C  gives  a  listing  of  the  subroutines  [8] 
which  were  used  to  solve  the  system  of  equations. 
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Equation  5  requires  the  surface  temperatures  as  known  data. 
These  temperatures  may  be  found  through  measurements  on  a  model  or  by  a 
complete  analysis  which  includes  all  modes  of  heat  transfer. 


2.2.  VIEW  FACTOR  EVALUATION 


The  solution  of  equation  5  requires  the  evaluation  of  the  view 
factors,  Fjs,  between  each  of  the  surfaces  in  the  enclosure.  Physically, 
the  view  factor  between  surface  i  and  surface  j  is  the  fraction  of  the 
total  radiant  energy  leaving  surface  i  which  strikes  surface  j  directly. 

The  view  factor  is  a  function  of  the  orientation  of  two  surfaces 
with  respect  to  each  other  and  the  distance  between  them.  The  view  factor 
between  two  finite  surfaces,  as  shown  in  Figure  2,  is  given  by  [1]: 


Fij  = 


1  //  CQ39i  C039j  dA.^ 


Ai 


irr 


ij 


(8) 


Solution  of  equation  8,  in  closed  form,  is 
number  of  surface  shapes  and  orientations  [1,2, 4, 5]. 
to  assemble  several  surfaces  with  known  view  factors 
surfaces  for  which  the  view  factor  may  be  determined 

El]. 


possible  for  a  limited 
It  is  often  possible 
to  form  irregular 
by  view  factor  algebra 


Equation  8  may  be  solved  by  direct  numerical  integration,  but 
this  requires  a  considerable  amount  of  computer  time.  Monte  Carlo  methods 
[2]  have  been  used  successfully  to  obtain  view  factors,  and  these  methods 
have  the  advantage  that  equation  8  does  not  need  to  be  solved  directly. 
Alternatively,  considerable  savings  in  computer  time  may  be  obtained  with 
direct  numerical  integration  of  the  view  factor  equation  by  first  applying 
Gauss'  law  to  equation  8  [6].  This  reduces  the  integration  from  a  double 
area  integral  to  a  double  contour  integral.  Application  of  Gauss'  law  to 
equation  8  yields: 

F,  ,  =  _L  /Mn(rn)  dSidSi  (9) 

1J  2irAi  1J  J  1 

where  r • ^  is  now  the  vector  between  points  on  the  contour  of  each  surface. 


Surface  j 


Figure  2.  Development  of  the  view  factor,  Fj-  ,  from  surface  "i"  to 
surface  "j".  (After  [1]). 
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View  factors  for  the  conical  tent  used  in  the  theoretical  portion 
of  this  study  can  be  found  by  knowing  the  view  factors  between  parallel 
concentric  disks,  and  applying  view  factor  algebra  to  obtain  the  view 
factors  between  actual  surfaces  (Appendix  B).  These  view  factors  should 
not  vary  appreciably  for  the  5-Man  Arctic  Tent  used  in  the  experimental 
portions  of  this  study.  A  closed  form  solution  to  equation  8  exists  for 
parallel  concentric  disks  [5]  and,  using  the  notation  of  Figure  3,  is  found 
to  be: 


Fjj  =  y2(z  -  /z"  -  Ux'y*) 

where 


(10) 


x  =  a/h 
y  =  h/b 

z  =  1  +  (1  +  x2)y2 


Evaluation  of  equation  9  by  numerical  integration  was  performed 
using  the  program  listed  in  Appendix  A.  Although  the  program  uses  a  simple 
integration  technique,  the  results  obtained  were  in  close  agreement  with 
those  obtained  using  the  exact  solution.  Table  1  lists  the  results 
obtained  for  different  step  sizes  and  different  disk  dimensions  along  with 
the  exact  solution  in  each  case.  The  error  incurred  by  using  a  relatively 
coarse  step  size  of  0.1  radians  was  less  than  7%  in  all  cases,  the  greatest 
error  occur ing  in  the  analysis  between  the  two  smallest  disks.  This  was  as 
expected  since,  as  the  disks  become  smaller,  or  closer  together,  a  fixed 
step  size  of  angle  around  the  disk  produced  an  arc  length  which  was  closer 
to  the  characteristic  dimension  of  the  problem.  Decreasing  the  step  size 
from  0.1  radians  to  0.01  radians  produced  only  marginal  improvements  while 
significantly  increasing  computation  time.  Increasing  the  step  size  to  0.5 
radians  caused  a  substantial  increase  in  the  error. 


Table  1 .  Comparison  of  results  of  view  factor  evaluation 
between  numerical  evaluation  of  equation  10  and 
the  exact  solution:  (a)  with  various  step  sizes; 
(b)  with  various  disk  dimensions.  Angles  are 


expressed  in  radians;  lengths 

are  expressed  in 

metres. 

(a)  3=0.254,  b=0 . 465 ,  h=0.27,  Fj. 

j  =0 . 21 2 

(exact) 

Step  Size 

A<J.  F 

Error 

% 

0.01  0.226 

+  6 

0.10  0.227 

+7 

0.50  0.243 

+  14 

8  - 


Surface  j 

I 

h 


V 


Surface  i 


Figure  3.  '  Configuration  of  the  concentric  disks  for  which  the 
view-factor  is  described  by  equation  10. 
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(b)  Step  Size,  A<j>=  0.10 


Disk 

a 

Dimensions 
b  h 

F  (comp) 

ij 

F  (exact) 

ij 

Error  % 

0.47 

0.93 

0.33 

0.217 

0.216 

0.5 

0.93 

1.17 

0.40 

0.513 

0.51 1 

0.4 

1.17 

1 .91 

0.45 

0.348 

0.348 

0.0 

1 .91 

1.91 

0.50 

0.774 

0.774 

0.0 

2.3  FABRIC  PROPERTIES 


The  thermal  radiative  properties  of  fabrics  have  been  found  to 
vary  considerably  between  different  materials  and  types  of  fabric 
construction  [3],  The  thermal  radiative  properties  of  several  fabrics 
which  are  commonly  used  in  tentage  are  listed  in  Table  2. 


Table  2.  Thermal  Radiative  Properties  Of  Selected  Tent 
Fabrics.  (Quoted  values  indicate  ranges  and 
[3]  should  be  consulted  for  actual  values  of 
specific  fabrics.) 


Material 

Mass/Area 

(kg/m2) 

Emissivity 

Reflectivity 

Transmissivity 

Nylon 

Cotton 

Cotton/Nylon 

Polyester 

Polypropylene 

0.05  -  0.11 
0.11  -  0.24 
0.19  -  0.30 
0.07  -  0.27 
0.07  -  0.15 

0.49  -  0.51 
0.60  -  0.82 
0.51  -  0.59 
0.54  -  0.85 
0.40  -  0.51 

0.33  -  0.41 
0.12  -  0.31 
0.33  -0.41 

0.11  -  0.38 

0.40  -  0.50 

0.06  -  0.18 
0.07  -  0.12 
0.07  -  0.08 
0.08  -  0.14 
0.07  -  0.18 

2.4  EXTERNAL  BACKGROUND  TEMPERATURE 


The  external  background  temperature  influences  the  internal 
thermal  radiant  exchange  through  the  transmissivity  of  the  tent  surfaces. 


10 


As  the  transmissivity  of  most  tent  fabrics  is  small,  and  since  multiple 
layers  of  fabric  are  often  used,  the  external  background  temperature  is 
usually  of  secondary  importance  in  determining  the  internal  radiant  heat 
fluxes.  It.  is,  however,  a  significant  factor  in  the  overall  radiant  heat 
transfer  from  the  tent.  For  typical,  single  layer  fabrics,  a  temperature 
difference  of  10  C  between  the  background  and  a  surface  inside  the  tent 
corresponds  to  an  additional  heat  transfer  of  approximately  2  W/m2  of 
radiant  energy  from  that  surface  to  the  external  background,  assuming  a 
transmissivity  of  0.15  for  the  tent  walls.  For  the  same  temperature 
difference  between  the  tent  fabric  and  the  external  background,  the  radiant 
energy  transfer  by  emission  from  the  fabric  will  be  approximately  7  to '13 
W/m2. 

t  ■ 

In  general,  the  background  temperature  will  be  some  combination 
of  the  temperatures  of  the  sky  and  ground,  as  well  as  their  respective 
radiative  properties  in  the  infrared  region  of  the  spectrum  and  the  view 
factors  from  the  external  tent  surfaces  to  the  sky  and  ground.  For  the  , 
purposes  of  this  study,  it  will  be  assumed  that  the  sky  dominates  the 
external  radiant  heat  flux  and  that  the  ground  has  the  same  temperature  as 
the  sky.  This  would  correspond  to  a  worst  case  analysis  for  the 
experiments  of  this  study  as  the  sky  was  found  to  be  slightly  colder  than 
the  ground.!  This  approximation  significantly  simplifies  the  analysis  as 
external  view  factors  for  the  tent  surfaces  to  the  surroundings  become  1. 
The  errors  incurred  by  using  this  approximation  are  small  for  this  study  as 
experimental  results  were  chosen  from  times  when  the  sky  temperature  was 
close  to  the  ambient  temperature.  This  approximation  would  result  in 
serious  errors  if  solar  heating  occurs  or  if  the  sky  temperature  varies 
substantially  from  the  ground  temperature.  External  view  factors  and  solar 
heating  will  be  studied  in  a  second  investigation  on  radiant  heat  transfer 
in  tents. 


Tyo  methods  of  predicting  the  sky  temperature  from  conventional 
meteorological  data  have  been  reported  [7].  The  simplest  of  these 
relationships  depends  only  upon  the  ambient  air  temperature  and  takes  the 
form: 


T$ky  ‘  °-°552  Tj-5 


(11  ) 


Another  relationship  includes  the  dew  point  temperature  in  ah  attempt  to 
account  for  moisture  in  the  air: 


TSky  =  Ta(  0.8  +  Tdp  ~  273  )1/4 

250 


(12) 
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In  the  experimental  portion  of  this  study,  a  pyranometer  was  used 
to  monitor  the  sky  temperature.  This  device  can  be  used  to  calculate  the 
sky  temperature  by  using  the  following  equation: 


oT 


•* 

sky 


3s  +  °Ts 


(13) 


where  it  has  been  assumed  that  the  sensor  and  background  emissivities  are 
both  1 . 


3.0  EXPERIMENT 


A  Canadian  Forces  5-Man  Arctic  Tent  was  set  up  in  an  open  field 
and  instrumented  with  thermocouples,  thermistors  and  heat  flow  disks.  Due 
to  the  geometric  symmetry  of  the  tent,  only  one  sector  of  the  tent  was 
instrumented.  Heat  was  supplied  by  two  forced  air  electric  heaters  with  a 
nominal  heat  output  of  1300  W  each. 

Wall  temperatures  were  measured  by  mounting  pairs  of 
thermocouples,  connected  in  parallel,  on  the  surface  of  the  tent  wall  as 
shown  in  Figure  4.  It  was  assumed  that  the  small  wire  size  (30  Gauge)  and 
the  pairing  of  the  thermocouples  would  improve  the  measurement  of  the  wall 
temperature.  Floor  temperatures  were  measured  by  differential 
thermocouples  at  each  of  the  heat  flow  disk  positions. 

Heat  flow  disks  were  placed  on  the  floor  of  the  tent  in  pairs  at 
four  positions.  The  disks  were  placed  so  as  to  give  an  equal  area 
weighting  of  the  heat  flux  to  the  floor.  For  each  pair  of  heat  flow  disks, 
one  was  covered  with  a  flat-black  paint  (e  =  0.95)  which  was  used  to 
measure  the  total  heat  transfer  rate  due  to  radiation  and  conduction  to  the 
floor.  The  other  heat  flow  disk  was  covered  with  a  layer  of  aluminum  foil 
(e  =  0.05),  attached  with  heat  sink  compound.  This  heat  flow  disk  was  used 
to  measure  the  conductive  component  of  heat  flow  to  the  floor.  The 
radiative  component  of  the  heat  flow  to  the  floor  was  calculated  from  the 
difference  between  these  two  measurements. 

Measurements  were  taken  at  night  to  eliminate  solar  heating, 
thereby  minimizing  the  effect  of  external  conditions  on  the  internal 
radiative  heat  fluxes.  Solar  heating  results  in  asymmetric  radiant  heat 
loss  transfer  to  the  tent  which  further  complicates  the  analysis.  The  sky 
temperature  was  calculated  from  a  measurement  with  a  pyranometer, 
referenced  to  the  air  temperature,  using  equation  13. 


Figure  M.  Configuration  of  the  thermocouples  used  for  measuring  the 
tent  wall  temperatures. 
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4.0  RESULTS 


The  output  voltage  of  the  thermocouples  and  the  heat  flow  disks 
was  small,  typically  0.1  mV  and  15  pV  respectively.  Signal  noise  with 
these  devices  was  found  to  be  approximately  ±5  pV .  The  signal  noise  was 
presumably  caused  by  electromagnetic  interference  picked  up  by  the  long 
lead  wires  from  the  tent  to  the  data  acquisition  system.  The  use  of  signal 
amplifiers  would  have  been  advantageous,  however,  none  was  available.  As  a 
result,  the  measurements  obtained  could  only  be  expected  to  give  an  order 
of  magnitude  confirmation  of  the  theoretical  predictions  of  the  radiant 
heat  transfer. 

Sky  temperatures,  as  determined  using  equation  11,  were  found  to 
underestimate  the  sky  temperature  evaluated  from  equation  13  from  between 

15  and  20  C,  or  by  10 %  of  the  absolute  temperature.  This  would  typically 
result  in  errors  in  the  radiant  heat  transfer  calculations  of  approximately 
4  W/m2  for  transmitted  radiant  energy  and  20  W/m2  for  emitted  radiant 
energy  by  the  fabric  surface. 

The  temperature  and  heat  flow  measurements  were  recorded  under 
uncontrollable  environmental  field  conditions.  Thus,  each  experimental  set 
of  results  was  unique.  Therefore,  the  complete  set  of  results,  nineteen  in 
all,  are  not  presented  here,  but  rather  typical  wall  and  floor  temperatures 
and  radiant  heat  fluxes  are  given,  as  shown  in  Figure  5.  Considering  the 
signal  noise  and  the  precision  of  the  reference  thermistors  (±0.2  C),  the 
wall  surface  temperature  readings  are  thought  to  be  accurate  to  within  ±0.5 
C.  Scatter  in  the  floor  surface  temperature  measurements  indicated  a 
variability  of  as  much  as  ±1  C  over  the  floor.  The  precision  of  the 
surface  temperature  measurement  corresponds  to  a  precision  in  the  computed 
radiant  heat  flow  between  15  and  5015  as  the  difference  between  the  wall 
temperatures  and  the  floor  temperature  varied  between  3  and  10  C  over  the 
series  of  experiments.  The  measured  radiant  heat  flux  from  the  floor  to 
the  tent  walls  was  found  to  be  between  20  and  60  W/m2  with  an  arithmetic 
average  of  all  measurements  of  approximately  40  W/m2. 

The  values  of  R^ji  f°r  each  experiment  was  calculated  using  the 
measured  values  of  surface  and  sky  temperatures  with  equation  5.  These 
values  of  radiant  heat  flux  were  then  used  with  equation  7,  and  the 
resulting  radiant  heat  flux  from  the  floor  computed.  The  computed  radiant 
heat  flux  to  the  floor  was  then  compared  with  the  measured  value  of  the 
relevant  experiment.  In  all  cases,  the  computed  radiant  heat  flux  from  the 
floor  was  less  than  the  measured  radiant  heat  flux  but  it  was  of  the  same 
order  of  magnitude.  Computed  radiant  heat  flows  from  the  floor  varied  from 

16  to  36  W/m2  with  a  numerical  average  of  25  W/m2.  Differences  between 
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Figure  5.  Typical  surface  temperatures  and  net  radiant  heat  fluxes 
for  the  5-Man  Arctic  Tent  when  heated  by  a  2000  Watt  heat 
source. 
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computed  and  measured  values  ranged  between  5  and  70$.  This  is  consistent 
with  the  temperature  measurement  precision  and  the  observed  signal  noise. 

The  temperatures  of  the  surfaces  of  the  tent  depend  upon  the 
ambient  temperature  and  the  rate  of  heating  of  the  tent.  Thus,  the  radiant 
heat  fluxes  found  in  this  study  are  specific  to  the  conditions  at  the  time 
of  the  experiments.  The  values  of  temperature  and  radiant  heat  flux  shown 
in  Figure  5  represent  typical  values  for  the  study  only  and  the  radiant 
heat  transfer  analysis  should  be  applied  for  cases  were  ambient  temperature 
or  rate  of  heating  differ.  It  should  also  be  noted  that  these  values  are 
specific  to  the  5-Man  Arctic  Tent  as  the  view  factors  for  tents  of 
different  shapes  and  sizes  will  be  different. 

Due  to  the  small  signal  of  the  heat  flow  disks,  the  observed 
signal  noise  and  possible  contamination  of  the  sensors  by  dust,  the 
measured  radiant  heat  flow  was  considered  to  be  less  accurate  than  that 
predicted  based  on  temperature  measurement. 

The  total  radiant  heat  loss  from  the  tent  to  the  external 
surroundings  was  found  to  be  approximately  600  to  750  W,  or  25$  of  the 
total  heat  loss  from  the  tent  in  this  case.  These  values  will  depend  upon 
the  tent  wall  construction,  the  amount  of  heat  supplied  to  the  tent  and  the 
tent  fabric  properties.  Total  radiant  heat  loss  from  the  tent  walls, 
internal  plus  external  radiant  heat  transfer,  was  found  to  be  approximately 
700  to  850  W. 


5.0  CONCLUSION 


Predictions  of  radiant  heat  transfer,  using  surface  temperature 
measurements  in  a  Canadian  Forces  5-Man  Arctic  Tent  with  a  model  tent  of 
comparable  dimensions,  agreed  with  measured  radiant  heat  transfer  rates  to 
the  floor  of  the  5-Man  Arctic  Tent  to  within  the  expected  precision  of  the 
experimental  method.  The  floor  was  found  to  be  a  major  contributor  to  the 
radiant  heat  transfer  in  the  tent.  The  magnitude  and  sense  of  the  radiant 
heat  transfer  from  the  floor  will  depend  upon  the  ground  temperature 
history  and  the  rate  of  heat  supplied  to  the  floor  of  the  tent. 

The  radiant  heat  loss  from  this  single  walled  tent  to  the 
external  surroundings  was  found  to  be  approximately  600  to  750  W  when  heat 
was  supplied  to  the  tent  at  a  rate  of  2600  W  or  approximately  25$  of  the 
total  heat  loss  from  the  tent.  Approximately  100  W  (net)  of  radiant 
thermal  energy  was  transferred  from  the  internal  surfaces  to  the  interior 
of  the  tent  which  was  subsequently  absorbed  by  the  walls  or  transmitted  to 
the  exterior  surroundings.  The  total  net  radiant  heat  transfer  from  all  of 
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the  surfaces,  both  external  and  internal,  was  found  to  be  700  to  850  W. 

This  quantity  -of  thermal  energy  must  be  supplied  to  the  walls  by  conductive 
transfer  from  the  heated  internal  air. 

Vi4w  factor  evaluation  by  numerical  integration  utilizing  Gauss' 
law  was  found  to  be  both  accurate  and  fast.  Computation  of  view  factors  is 
made  easier  by  using  interstitial  geometric  surfaces  which  are  easily 
described  mathematically  in  computer  code  and  applying  view  factor  algebra 
to  obtain  the  view  factor  between  desired  surfaces  and  the  interstitial 
surfaces.  j 

Sky  temperature  prediction  based  on  ambient  temperature  alone  was 
found  to  be  accurate  to  approximately  10$  for  overcast  skies  or  warm,  clear 
skies.  This  results  in  an  error  in  the  radiant  heat  transfer  computations 
of  approximately  20  W/m2  for  net  emitted  radiant  energy  and  4  W/m2  for 
radiant  energy  transmitted  by  a  single  fabric  layer. 
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2: 

3:  C 
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5  C 
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12:  C 
13:  C 
14:  C 
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16:  C 
17:  C 
18:  C 
19:  C 
20:  C 
21  :  C 
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23:  C 
24*  C 
25*  C 
26:  C 
27:  C 
28:  C 
29.  C 
30:  C 
31  :  C 
32*  C 
33:  C 
34:  C 
35:  C 
36:  C 
37:  C 
38:  C 
39 :  C 
40:  C 
41  :  C 
42:  C 
•43:  C 
44*  C 
45* 

46* 

47* 

48*  C 
49* 

50* 

51*  10 

52*  C 
53* 


A.  Computer  Program  For  Evaluation  Of  View  Factors  Between 
Concentric  Disks.  *  • 

PROGRAM  CVIEW 


THtS  PROGRAM  IS  INTENDED  FOR  DETERMINING  THE 
VlEWFACTOR  BETWEEN  TWO  SURFACES  FOR  USE  IN  RADIATIVE 
HEAT  TRANSFER  PROBLEMS.  THE  PROGRAM  MAKES  USE  OF  THE 
CONTOUR  INTEGRAL  FORMULATION  INSTEAD  OF  THE 
CONVENTIONAL  AREA  INTEGRAL. 

THE  INFORMATION  ON  THIS  TECHNIQUE  MAY  BE  FOUND  IN: 

Sparrow ;A  new  and  simpler  formulation  for  radiative 
f  _  ANGLE  FACTORS; J -OF  HEAT  TRANS. ;C85 ;MAY 1963 • 

IT  USES  A  TRAPEZIODAL  RULE  TO  DO  THE  INTEGRATION  THE 
LIMITS  OF  THE  CONTOURS  AROUND  THE  TWO  SURFACES.  THE 
.AREA  OF  SURFACE  "1"  AND  THE  INCREMENTAL  STEP  SIZES  ARE 
READ  IN  AS  DATA.  THE  DIRECTION  OF  THE  INTEGRATION 
AROUND  THE  CONTOURS  IS  SUCH  THAT  IF  AN  OBSERVER  WERE 
TO  WALK  AROUND  THE  CONTOUR  IN  THE  DIRECTION  OF  THE 
INTEGRATION  WITH  HIS  HEAD  POINTING  IN  THE  DIRECTION 
OF  THE  SURFACE  NORMAL  AT  THE  CONTOUR  THEN  THE  AREA 
BOUNDED  BY  THAT  CONTOUR  WOULD  ALWAYS  LIE  TO  THE 
LEFT  OF  THE  OBSERVER. 

THIS  UPDATE  PROGRAM  IS  FOR  USE  WITH  CIRCULAR  CONTOURS 
LYING  IN  THE  R-THETA  PLANE. 

TH?  CORRECT  SOLUTION  TO  THIS  PROBLEM  MAY  BE  FOUND  IN: 
OZISIK;  BASIC  HEAT  TRANSFER;  MCGRAW-HILL 

AREA1  =  THE  AREA  OF  THE  FIRST  SURFACE 
LIMIT ()  =  THE  LIMITS  OF  THE  CONTOURS  AROUND  THE 
TWO  SURFACES 

R , THETA, Z  =  THE  POLAR  COORDINATES  OF  THE  CONTOURS 
F( )'  =  THE  VALUE- OF  THE  INTEGRAND  AT  THE  PRESCRIBED 

j  COORDINATES 

SUM  =  THE  RUNNING  SUM  OF  THE  INTEGRATION 
HTRAP  .  =  THE  VALUE  OF  THE  MULTIPLICATION  OF  THE 
i  SEP  SIZES. 

FTRAP  =  THE  VALUE  OF  THE  INTEGRAL  FOR  SPECIFIC 
COORDINATES 

ITRAP  =  THE  VALUE  OF  THE  INTEGRAL  OVER  A  PORTION 
OF  THE  CONTOUR. 

F 12  =  THE  VALUE  OF  THE  VlEWFACTOR  BETWEEN 

SURFACES  ONE  AND  TWO. 

REAL  LIMIT (12) .THETASTEP . THETA 1STEP .THETA2STEP , PI 
REAL  R 1 , R2 , Z 1 , Z2 , FTRAP , F ( 4  V, HTRAP , ITRAP , SUM , AREA  1 
REAL  THETA  IP ,THETA2P 

DO  10  1=1’,  12 

LIMIT (I )=0 . 0  ■  • 

CONTINUE 

PI=3- 141592654 
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54* 

55*  C 
56* 

57* 

58*  C 
59* 

60* 

61*  C 
62*  C 
63*  C 
64*  C 
65*  C 
66*  C 
67* 

68* 

69*  C 
70* 

71*  C 

72* 

•73* 

74*  C 
75* 

76* 

77*  C 
78* 

79* 

80* 

81* 

82* 

83*  C 
84* 

85* 

86*  20 
87*  C 
88* 

89* 

90*  C 
91*  30 

92*  40 

93: 

94:  C 
95: 

96:  190 
97: 

98:  C 
99: 

100: 


SUM  =  00 


LIMIT(2)=2*PI 

LIMIT(3)=2*PI 

READ( 105)*)  R1 ,R2,Z1 , Z2,THETASTEP 
AREA  1 =PI*R 1**2 


THIS  PORTION  OF  THE  PROGRAM  IS  SET  UP  TO  DETERMINE 
THE  VIEWFACTOR  BETWEEN  TWO  CONCENTRIC  HORIZONTAL 
DISKS. 

THETA  1 STEP=THETASTEP 
THETA2STEP=-THETASTEP 

HTRAP=THETA 1STEP*THETA2STEP 

DO  40  THETA  1 =LIM IT ( 1 ) , LIMIT (21 , THETA  1  STEP 
THET A  1 P  =THET A 1 +THETA 1 STEP 

DO  30  THETA2=LIMIT(3)  LIMIT (4)  THETA2STEP 
THETA2P =THETA2+THETA2STEP 

F ( 1 ) =FUN (R 1 . THETA  1 , Z 1 , R2  THETA2 , Z2 ) 

F(2)=FUN (R  1 , THETA  IP ; Z  1 , R2  THETA2 , Z2) 
F(3)=FUN(R1 , THETA  1 ,Z1 , R2 , THETA2P , Z2 ) 

F(4 )=FUN (R  1 : THETA  IP : Z 1 , R2 ,THETA2P  Z2) 
FTRAP=0.0 


DO  20  1=1,4 

FTRAP=FTRAP+F ( I ) 

CONTINUE 

ITRAP=HTRAP*FTRAP/4 .0 
SUM=SUM+ITRAP 

CONTINUE 

CONTINUE 

F12=SUM/( AREA  1*3- 14159*2) 

W R ITE ( 6  190)  F 12 , AREA  1 

FORMAT (  0  /THE  VIEWFACTOR  FROM  SURFACE  1  TO  2  IS: 
+  F 1 0  •  7 ,  /  /  THE  AREA  OF  SURFACE  1  IS:  /F10-3) 

RETURN 

END 
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1 :  C 

2:  FUNCTION  FUN (X 1 , Y 1 , Z 1 , X2 , Y2 , Z2 ) 

3:  C 

4 :  c*****  THIS  FUNCTION  EVALUATES  THE  INTEGRAND  FOR  THE 
5  CONTOUR  INTEGRAL  AT  THE  PRESCRIBED  POINTS  ON 

6;  c****»  THE  CONTOUR. 

7:  C  : 

8*  C 

9*  REAL  R,X1,X2,Y1,Y2,Z1,Z2 

io*  c  !  v' 

11*  C  NOW,  X,  Y.  Z  REFER  TO  R, THETA ,  Z  RESPECTIVELY 
12*  C 

13*  COSY  =  C0S(Y1-Y2) 

14*  C 

15*  R=SQRT(X 1 **2+X2**2  2*X 1 *X2*C0SY+(Z2  Z1)**2) 

16*  FUN  =  L0G(R)*X1*X2*C0SY 

17:  C 

1 8 :  RETURN 

19:  END 


APPENDIX  B 


EXAMPLES  OF  VIEW  FACTOR  ALGERBRA 


It  is  assumed  that  the  view  factors  between  surfaces  1,2  and  3  of 
Figure  B-1  are  known  by  evaluating  the  view  factor  equation  for  these 
concentric  disks  (so  chosen  for  the  ease  of  evaluation  of  the  view  factor 
equation).  From  conservation  of  energy,  the  sum  of  all  view  factors  from  a 
surface  must  equal  1 .  The  view  factor  from  surface  3  to  surface  5  can  then 
be  found  from: 


F35  +  F32  =  1  (B1 ) 

Similarly: 

F34  +  F35  +  F31  =  1  (B2) 


Alternatively: 


F34  "  f32  x  F24 


(B3) 


Finding  one  view  factor  between  two  surfaces  and  knowing  the 
surface  areas  of  the  two  surfaces  provides  sufficient  information  to 
determine  the  remaining  unknown  view  factor: 


A5  x  F53  “  A3  x  F35 


(B4) 


as  this  relationship  is  required  for  conservation  of  radiant  thermal 
energy. 


IMAGINARY  SURFACE 


Figure  B1 .  Example  configuration  for  view-factor  algebra. 
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APPENDIX  C 


COMPUTER  PROGRAM  FOR  SOLUTION  OF  THE  RADIATIVE 
HEAT  TRANSFER  IN  TENTS 


Page 

1 .  RADQ . 24 

2.  DOOMS.  .  .  . .  26 

3.  SOLVE .  29 

4.  IMPRV . . .  .  30 
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1 

1  , 
1 
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Appendix  C.  Computer  Program  For  Solution  Of  The  Radiative  Heat 
Transfer  In  Tents. 


1  C 

2 

3  C 

4  C 

5  C 

6  C 

7  C 

8  C 

9  C 
10 

1  1 
12 

13 

14 

15 

16  C 

17  C 

18 

19 

20 
21 

22  C 

23  C 

24  C 

25 

26  C 

27 

28  C 

29 

30 

31 

32 

33 

34  C 

35 

36 

37  C 

38  10 

39  C 

40  20 

41  C 

42  C 

43  C 

44 

45  C 

46  C 

47  c 

48 

49  C 

50 

51 

52 

53  C 

54  C 


PROGRAM  RADQ 


THIS  PROGRAM  IS  INTENDED  TO  SOLVE  THE  SYSTEM  OF 
EQUATIONS:  M.R  =  T,  TO  DETERMINE  THE  RADIOSITIES 
FROM  THE  SURFACES  OF  AN  ENCLOSURE 

REAL  M(20 , 20)  , RLI (20 ) , T (20 >  ,TAMB,RHO(20)  ,EPS(20'> 
REAL  SIGMA, DQ,D(20) ,MM(20,20) ,Z(20) ,RES(20) 

REAL  FIJ(20,20) ,QME AS , TSKY , QSKY , DELT A ( 2 0 , 20 ) 

REAL  ROUT (.20)  ,RIN(20)  ,RRI(20')  ,TT(20)  ,T0UT(20) 
REAL  QNET(20) 

INTEGER  NPIV(20) , IND , L,N ,EXPT , EQ 


READC105,*)  L,N 

READ( 1 05, *)  ( (FIJ(I , J) ,J=1 ,N) ,1=1 ,N) 

READO05,*)  ( (EPS( I)  ,RH0(I ) )  ,  1=  1  ,N) 

SIGMA  =  5-6696E-8 

SET  UP  THE  STIFFNESS  MATRIX  OF  VIEWFACTORS  AND  RHO'S 
DO  20  1=1, N 
DO  10  J=  1 ,  N 

IF(I.EQ. J)THEN  .  . 

•  DELTA( I , J)  =1.0 
ELSE 

*  DELTA(I, J)  =0.0 
|  END  IF 

;  M( I , J)  =  DELTA ( I , J)  -  RHO(I)*FIJ(I, J) 

■  MM(I, J)  =  M( I , J) 

CONTINUE 

CONTINUE 

DO  THE  LU  DECOMPOSITION  OF  THE  STIFFNESS  MATRIX 
CALL  DC0MS(20,N,M,NPIV,D, IND) 

READ- IN  THE  TEMPERATURE  DATA  (REFERENCED  TO ' AMBIENT ) 
DO  66  K=  1  ,L 

RE AD [ 1 05 , * )  QMEAS , EXPT , T AMB , QSKY 
READt 105,*)  (T(I) ,1=1 ,N) 

TSKY  =  QSKY  +  SIGMA* (TAMB+273- 16)**4 

CONVERT  TEMPERATURES  TO  RADIANT  ENERGIES 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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1 

1 
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1 

1 

1 

1 
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1 

1 

1 

1 

1 
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55  C 

56 

57  c 

58 

59  C 

60  C 

61  C 

62 

63 

64 

65 

66 

67  C 

68  C 

69  C 

70 

71  C 

72  30 

73  C 

74  C 

75  C 

76 

77 

78  C 

79 

80 
81 
82 

83 

84 

85 

86  C 

87 

88  C 

89  C 

90  C 

91 

92  C 

93 

94  C 

95 

96  C 

97  40 

98  C 

99  C 

100  C 

101  C 

102 

103 

104  C 

105 

106  C 

107 

108  C 


DO  30  1=1 ,N 

T  ( I )  =  SIGMA*  (T  (D+TAMB+273 . 1 6  )**4 

EVALUATE  THE  AMBIENT  BACKGROUND  RADIANT  ENERGIES 

IF  (I.EQ.1)  THEN 
TOUT ( I )  =  T(I) 

ELSE 

TOUT (I)  =  TSKY 
END  IF 

CONSTRUCT  THE  RIGHT  HAND  SIDE  OF  THE  SYSTEM  OF  EQ. 

TT( I )  =  EPS(I)*T(I)  +  ( 1 -EPS( I )-RHO( I ) ) *TOUT ( I ) 

CONTINUE 

SOLVE  FOR  RLI ,  THE  INWARD  FLUX  FROM  THE  SURFACES 

CALL  SOLVE (20,N,M,NPIV,TT,RLI) 

CALL  IMPRV(20,N,MM,M,NPIV,TT,RLI,Z,RES1 

PRINT  *,'  ' 

PRINT  * ,  'EXPERIMENT  NUMBER  :  ',EXPT 
PRINT  */  ' 

PRINT  *,'  ' 

PRINT  * /SURFACE  NET  HEAT  NET  HEAT  NET  HEAT ' 

PRINT  */  LOSS  AMBIENT  LOSS  INTERNAL  LOSS  TOTAL' 

PRINT  */  ' 

DO  50  1=1, N 

SUM  UP  THE  INTERNAL , INBOUND  ENERGIES  TO  EACH  SURFACE 
RRI(I)  =0.0 
DO  40  J=1 ,N 


RRI(I)  =  RRI(I)  +  FIJ( I , J)*RLI( J ) 


CONTINUE 

COMPUTE  THE  NET  HEAT  LOSSES  INTERNALLY,  EXTERNALLY 
AND  TOTAL  FOR  EACH  SURFACE 

ROUT ( I )  =  EPS(I1*T(I)  -  (1 .-RHO(I))*TOUT(I)  + 
(1.-EPS(I)-RH0(I)}*RRI(I> 


RIN(I)  =  RLKII  -  RRI(I) 
QNET(I)  =  ROUT ( I )  +  RIN(I) 
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1 

109 

C 

PRINT  OUT  THE  RESULTS 

1 

1 10 

C 

1 

1 1  1 

C 

1 

1 12 

WRITE(106, 100)  I , ROUT ( I ) ,RIN(I) , QNET ( I ) 

1 

113 

100 

FORMAT ( '  ',3X,I2,3(F15.2)) 

1 

114 

C 

t. 

1 

1 15 

50 

CONTINUE 

1  116  C 

1  117  C 

1  118  C 

1  119  C 

1  120 

1  121 
1  122 

1  123 

1  124 

1  125 

1  126 

1  127 

1  128 

1  129 

1  130 

1  131 

1  132  C 

1  .133  C 

1  134  60  CONTINUE 

1  135  STOP 

1  136  END 

2  1  c  : 

2  2  C  ! 

2  3  SUBROUTINE  DCOMS(NDIM, N , A, NPIV , D, IND1 

2  4-  C 

2  5  C 

2  ,  6  C  THIS  SUBROUTINE,  DECOMPOSITION  WITH  SCALED  PARTIAL  PIVOTING, 

2  7  C  DOES  GAUSSIAN  ELIMINATION  OR,  EQUIVALENTLY,  A  TRIANGULAR  (LU) 

2  8  C  FACTORIZATION  OF  THE  N*N  MATRIX  STORED  IN  THE  ARRAY  "A".  AT 

2  9  C  COMPLETETION ,  THE  "A"  WILL  CONTAIN  THE  LOWER  TRIANGULAR  MATRIX- 

2  10  C  OF  MULTIPLIERS  USED  IN  THE  ELIMINATION  AS  WELL  AS  THE  UPPER 

2  11C  TRIANGULAR  MATRIX  "U",  THE  RESULT  OF  THE  ELIMINATION. 

2  12  C  : 

2  13  C  THE  MATRIX  IS  ASSUMED  TO  BE  SINGULAR  IF  EITHER  SOME  ROW  IS 

2  14  C  ZERO  INITIALLY,  OR,  SOME  "SCALED"  PIVOT  DURING  THE  ELIMINATION 

2  15  C  IS  SMALLER  THAN  UNIT  ROUND-OFF.  IF  THE  FORMER  HOLDS,  THE 

2  16  C  DECOMPOSITION  DOES  NOT  COMMENCE.  IN  THE  LATTER  CASE,  DCOMS 

2  17  C  WILL  COMPLETE  THE  DECOMPOSITION  BUT  THE  RESULTING  UPPER 

2  18  C  TRIANGULAR  MATRIX  WILL  BE  SINGULAR. 

2  19  C  I 

2  20  C  THIS  ROUTINE  CAN  BE  USED  IN  CONJUNTION  WITH  THE  ROUTINE 

2  21  C  SOLVE  TO  FIND  THE  SOLUTION  TO  A  SYSTEM  OF  LINEAR  EQUATIONS. 

2  22  C  THE  .CURRENT  MAXIMUM  SIZE  OF  THE  SYSTEM  IS  50*50 

2  23  c  ; 

2  24  C  TAKEN  FROM:  STUNT  MANUAL;  DEPT  OF  COMPUTER  SCIENCE ;U  OF  T 

2  25  c  ; 

2  26  C  ' 


PRINT  OUT  COMPARISON  BETWEEN  COMPUTED  AND  MEASURED 
RADIANT  HEAT  TRANSFER  AT  THE  FLOOR 

PRINT  ' 


PRINT  *, 'COMPUTED  RADIATIVE 
PRINT  * , 'MEASURED  RADIATIVE 
PRINT 

DQ  =.QNET(1)  -  QMEAS 
PRINT  *,' 

EQ  =  (  DQ/QMEAS  ) *  1 00 
PRINT  *,  ' 

PRINT  *,  '  ' 

PRINT  *,'  ' 

PRINT  *,  '  ' 

print  *,'  : 


FLOOR  HEAT  LOSS  =  ' ,QNET( 1 ) 
FLOOR  HEAT  LOSS  =  ' , QMEAS 


DIFFERENCE  =  '  ,DQ 

ERROR  =  '  ,EQ, '%' 
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2 

2?  C 

CALLING 

SEQUENCE:  CALL  DCOMS(NDIM,N, A, NPIV, D, IND) 

2 

28  C 

2 

29  C 

PARAMETERS: 

2 

30  C 

2 

31  C 

NDIM 

-  AN  INTEGER  INDICATING  THE  NUMBER  OF  ROWS 

2 

32  C 

IN  THE  ARRAY  "A"  AS  DECLARED  IN  THE 

2 

33  C 

CALLING  PROGRAM. 

2 

34  C 

2 

35  C 

N 

-  AN  INTEGER  CONSTANT  INDICATING  THE  SIZE  OF 

2 

36  C 

THE  SYSTEM  TO  BE  SOLVED 

2 

37  C 

2 

38  C 

A 

-  A  REAL  2-DIMENSIONAL  ARRAY,  OF  SIZE  NDIM*N 

2 

39  C 

HOLDING  THE  MATRIX  TO  BE  DECOMPOSED.  ON 

2 

40  C 

RETURN,  THE  CONTENTS  OF  "A"  ARE  REPLACED 

2 

41  C 

BY  THE  LU  FACTORIZATION. 

2 

42  C 

2 

43  G 

NPIV 

-  AN  INTEGER  VECTOR,  OF  SIZE  N,  WHICH  IS 

2 

44  C 

UNINITIALIZED  AT  THE  TIME  OF  CALLING.  THIS 

2 

45  C 

ARRAY  WILL  RECORD  THE  REARRANGING  OF  THE 

2 

46  C 

2 

47  C 

D 

-  A  REAL  VECTOR,  OF  SIZE  N,  THAT  IS  USED  AS 

2 

48  C 

A  WORKSPACE  FOR  THE  SCALING  OPERATION.  THIS 

2 

49  C 

ARRAY  IS  DECLARED  BUT  NOT  INITIALIZED  IN 

2 

50  C 

THE  CALLING  PROGRAM. 

2 

51  C 

2 

52  C 

IND 

-  AN  INTEGER  INDICATING  IF  "A"  IS  SINGULAR  OR  NOT 

2 

53  C 

IND  =  -1  ;  "SCALED"  PIVOT  <  UNIT  ROUND-OFF 

2 

54  C 

=  0  ;  "A"  IS  NONSINGULAR 

2 

55  C 

THIS  PARAMETER  IS  NOT  INITIALIZED  AT  THE  TIME 

2 

56  C 

OF  CALLING. 

2 

57  C 

2 

58  C 

2 

59 

DIMENSION  A(NDIM,NDIM) ,D(NDIM) ,NPIV(NDIM) 

2 

60 

IND  = 

0 

2 

61  C 

2 

62  C 

CHECH  FOR  A  SYSTEM  OF  ONLY  ONE  UNKNOWN 

2 

63  C 

2 

64 

IF  (N 

.EQ.1)  RETURN 

2 

65  C 

2 

66  C 

INITIALIZE  PIVOT  AND  D  VECTORS 

2 

67  C 

2 

68 

DO  20 

1=1  ,  N 

2 

69 

NPIV(I)  =  I 

2 

70 

D(I)=0.E0 

2 

71  C 

2 

72  C 

THE  LARGEST  ABSOLUTE  VALUES  IN  EACH  ROW  ARE  RECORDED  IN  "D" 

2 

73  C 

2 

74 

DO 

10  J=1 ,N 

2 

75 

D ( I )  =  AM AX  1 (D( I ) , ABS (A(I,J))) 

2 

76  10 

CONTINUE 

2 

77  C 

2 

78 

IF 

CD ( I )  .EQ.O.EO)  D(I)  =  1 .EO 

2 

79  20 

CONTINUE 

2 

80  C 
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2 

81 

C 

MAIN  LOOP  FOR  GAUSS  ELIMINATION 

2 

82 

C 

2 

83 

NM1  =  N-1 

2 

84 

C 

2 

85 

DO  :80  1  =  1  ,  NM1 

2 

86 

C 

i 

2 

87 

C 

DETERMINE  THE  LARGEST  "SCALED  PIVOT,  IE, 

2 

88 

C 

MAX  |A(J, I)/D(J) ! ,  I<=J<=N 

2 

89 

C 

‘ 

2 

90 

jCOLMAX  =  0.0 

2 

91 

C 

i 

2 

92 

DO  30  J=I,N 

2 

93 

IP  =  NPIV(J) 

2 

94 

HOLD  =  ABS(A(IP, I))/D(IP) 

2 

95 

IF  ( HOLD . LE . COLM AX )  GOTO  30 

2 

96 

COLMAX  =  HOLD 

2 

97 

NROW  =  J 

2 

98 

30 

CONTINUE 

2 

99 

C 

2 

100 

C 

TEST  FOR  SINGULARITY.  THE  MATRIX  IS  ASSUMED  TO  BE  SINGULAR 

2 

101 

C 

IF 

COLMAX  (THE  ABS.  VALUE  OF  THE  SCALED  PIVOT)  IS  EQUIVALENT 

2 

102 

C 

TO 

ZERO,  IE,  1.0  +  COLMAX  =  1.0 

2 

103 

C 

2 

104 

C 

IF 

THIS  IS  TRUE  THEN  THE  ROUTINE  PROCEEDS  ON  TO  THE  (1+1  )TH 

2 

105 

C 

STAGE (OF  THE  ELIMINATION 

2 

106 

c 

2 

107 

, : IF  ((1 .O+COLMAX)  .NE.  1.0)  GOTO  40 

2 

108 

IND  =  -1 

2 

109 

i  GOTO  80 

2 

1  10 

c 

* 

2 

1 1 1 

c 

IF 

AN  INTERCHANGE  IS  NECESSARY,  ALTER  THE  PIVOT  VECTOR  "NPIV" 

2 

M2 

c 

2 

M3 

40 

. IPIVOT  =  NPIV (NROW) 

2  . 

114 

'IF  (NROW  .EQ.  I)  GOTO  50 

2 

:  1  15 

NPIV (NROW)  =  NPIV(I) 

2 

1 16 

NPIV(I)  =  IPIVOT 

2 

1  17 

C 

2 

1 18 

C 

THE  MULTIPLIERS  FOR  THE  COMPUTATION  OF  THE  REMAINING  ROWS  ARE 

2 

119 

C 

DETERMINED  AND  ELINIMATION  IS  PERFORMED.  THE  VALUE  OF  EACH 

2 

120 

C 

MULTIPLIER  IS  STORED  IN  THE  POSITION  OF  THE  ELIMINATED 

2 

121 

c 

ELEMENT . 

2 

122 

c 

: 

2 

123 

50 

i  ip i  =  i+i 

2 

124 

'  DO  70  J=IP1 ,N 

2 

125 

JPIVOT  =  NPIV ( J) 

2 

126 

• 

;  AMULT  =  A( JPIVOT , I)/A( IPIVOT , I) 

2 

127 

|  A( JPIVOT, I)  =  AMULT 

2 

128 

C 

! 

2 

129 

‘  DO  60  K  =  IP  1 , N 

2 

130 

A( JPIVOT , K)=A( JPIVOT , K)-AMULT*A( IPIVOT , K) 

2 

131 

60 

CONTINUE 

2 

132 

70 

,  CONTINUE 

2 

133 

80 

'CONTINUE 

2 

134 

C 
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2  135  TEST  =  1.  +  ABS(A(NPIV(N) ,N)) 

2  136  IF  (TEST  .EQ.  1.)  THEN 

2  137  IND  =  -1 

2  138  END  IF 

2  139  C 

2  1 40  RETURN 

2  141  END 


3 

1 

C 

3 

2 

SUBROUTINE  SOLVE (NDIM , N , LU , NPIV , B , X ) 

3 

3 

C 

3 

4 

c 

3 

5 

c 

THIS  SUBROUTINE  PERFORMS  THE  FORWARD  AND  BACKWARD  SUBSTITUTIONS 

3 

6 

c 

STEPS  IN  THE 

SOLUTION  OF  A  SYSTEM  OF  LINEAR  EQUATIONS  AX=B. 

3 

7 

c 

IT  ASSUMES  THAT 

THE  TRIANGULAR  (OR  LU)  FACTORIZATION  OF  A  HAS 

3 

8 

c 

ALREADY  BEEN 

COMPUTED  BY,  SAY,  THE  ROUTINE  DCOMP  OR  DCOMS.  IF 

3 

9 

c 

EITHER  ROUTINE 

INDICATES  THAT  "A"  IS  SINGULAR,  THEN  THE  USE  OF 

3 

10 

c 

SOLVE  MAY  PRODUCE  AN  OVERFLOW  INDICATION. 

3 

1 1 

c 

3 

12 

c 

3 

13 

c 

REPRODUCED 

FROM:  STUNT  MANUAL,  DEPT  OF  COMPUTER  SCIENCE,  U  OF  T 

3 

14 

c 

3 

15 

c 

3 

16 

c 

CALLING  ; 

SEQUENCE:  CALL  SOLVE (NDIM, N , LU , NPIV , B , X) 

3 

17 

c 

3 

18 

c 

PARAMETERS: 

3 

19 

c 

3 

20 

c 

NDIM 

- 

AN  INTEGER  INDICATING  THE  NUMBER  OF  ROWS  IN  THE 

3 

21 

c 

ARRAY  "A"  AS  DECLARED  IN  THE  CALLING  PROGRAM 

3 

22 

c 

3 

23 

c 

N 

- 

AN  INTEGER  CONSTANT  INDICATING  THE  SIZE  OF  THE 

3 

24 

c 

SYSTEM  TO  BE  SOLVED. 

3 

25 

c 

3 

26 

c 

LU 

_ 

A  REAL  2-D  ARRAY  OF  SIZE  NDIM*N  CONTAINING  THE 

3 

27 

c 

LU  DECOMPOSITON  OF  A.  THIS  ARRAY  IS  NOT  ALTERED 

3 

28 

c 

BY  SOLVE. 

3 

29 

c 

3 

30 

c 

NPIV 

- 

AN  INTEGER  VECTOR,  OF  DIMENSION  N,  HOLDING  THE 

3 

31 

c 

PIVOT  INFORMATION  FOR  THE  ELIMINATION  STEP. 

3 

32 

c 

3 

33 

c 

X 

- 

A  REAL  VECTOR,  OF  SIZE  N,  THAT  IS  DECLARED  BUT 

3 

34 

c 

NOT  INITIALIZED  BY  THE  CALLING  PROGRAM.  ON 

3 

35 

c 

RETURN,  THIS  ARRAY  CONTAINS  THE  COMPUTED 

3 

36 

c 

SOLUTION  OF  THE  SYSTEM. 

3 

37 

c 

3 

38  C 

B 

- 

A  REAL  VECTOR,  OF  SIZE  N,  HOLDING  THE  RIGHT  HAND 

3 

39 

c 

SIDE  OF  THE  ORIGINAL  SYSTEM  TO  BE  SOLVED.  THE 

3 

40 

c 

CONTENTS  OF  THIS  VECTOR  ARE  UNALTERED  BY  SOLVE. 

3 

41 

c 

3 

42 

c 
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3 

43 

DIMENSION  B(NDIM) ,X(NDIM) ,NPIV(NDIM) 

3 

44 

REAL  LU(NDIM,NDIM) 

3 

45 

C 

3 

46 

C 

CHECK  FOR  SYSTEM  OF  ONLY  ONE  UNKNOWN 

3 

47 

C 

3 

48 

IF  ■(  N .  GT .  1 )  GOTO  10 

3 

49 

X(1)  =  B( 1 )/LU( 1,1) 

3 

50 

RETURN 

3 

51 

C 

3 

52 

c 

FORWARD  ELIMINATION  ON,  "B".  THE  RESULT 

3 

53 

C 

3 

54 

10 

KPIVOT  =  NPIV(I) 

3 

55 

X( 1 )=B(KPIVOT) 

3 

56 

DO  !30  K=2,N 

3 

57 

KPIVOT  =  NPIV(K) 

3 

58 

KM1  =  K-1 

3 

59 

!SUM  =  B(KPIVOT) 

3 

60 

DO  20  J=1 ,KM1 

3 

61 

i  •  SUM  =  SUM  -  LU (KPIVOT , J)*X( J) 

3 

62 

20 

CONTINUE 

3 

63 

;X(K)  =  SUM 

3 

64 

30 

CONTINUE 

3 

65 

C 

i  ■* 

3 

66 

c  : 

BACK  SUBSTITUTION  BEGINS 

3 

'  67 

c 

3 

68 

XQO  =  X  ( N )  /LU  (KPIVOT ,  N ) 

3 

69 

K  i  N 

3 

70 

DO  !50  1=2, N 

3 

71 

KP 1  -  K 

3 

72 

K  =  K-1 , 

3 

73 

•KPIVOT  =  NPIV(K) 

3 

74 

|SUM  =  X(K) 

3 

75 

iDO  40  J  =  KP1.N 

3 

76 

SUM  =  SUM  -  LU (KPIVOT, J)*X( 

3 

•  77 

40 

CONTINUE 

3 

78 

X(K)  =  SUM/LU (KP IVOT , K ) 

3 

•79 

50 

CONTINUE 

3 

80 

C 

3 

81 

RETURN 

3 

82 

END 

4  2  C 

4  3  SUBROUTINE  IMPRV(NDIM,N, A,LU,NPIV,B,X,Z,R1 

4  4  C 

4  5  C  ; 

4  6  C  GIVEN  AN  APPROXIMATE  SOLUTION  FOR  A  LINEAR  SYSTEM  OF  EQUATIONS 

4  7  C  THIS  SUBROUTINE  CARRIES  OUT  ONE  ITERATION  OF  THE  ITERATIVE 

4  8  C  IMPROVEMENT  PROCESS  FOR  COMPUTING  A  BETTER  APPROXIMATE  SOLUTION. 

4  9C  IT  IS  ASSUMED  THAT  BOTH  THE  ORIGINAL  MATRIX  AND  ITS  LU 
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4  10  C  DECOMPOSITION  ARE  AVAILABLE.  THIS  SUBROUTINE  CAN  BE  USED 

4  1 1  C  IN  CONJUNCTION  WITH  THE  SUBROUTINES  DCOMP,  DCOMS  AND  SOLVE 


4 

12  C 

TO  FIND  THE  SOLUTION  OF  A  LINEAR  SYSTEM  OF  EQUATIONS. 

4 

13  C 

4 

14  C 

TAKEN  FROM 

:  STUNT  MANUAL,  DEPT  OF  COMPUTER  SCIENCE,  U  OF  T 

4 

15  C 

4 

16  C 

CALLING  SEQUENCE:  CALL  IMPRV(NDIM,N , A, LU , NPIV ,B , X, Z , R) 

4 

17  C 

4 

18  C 

PARAMETERS 

: 

4 

19  C 

4 

20  C 

NDIM 

-  AN  INTEGER  CONSTANT  INDICATING  THE  NUMBER  OF 

4 

21  C 

ROWS  IN  THE  ARRAYS,  "A"  AND  "LU"  AS  DECLARED 

4 

22  C 

IN  THE  CALLING  PROGRAM 

4 

23  C 

4 

24  C 

N 

-  AN  INTEGER  CONSTANT  INDICATING  THE  NUMBER  OF 

4 

25  C 

UNKNOWNS  IN  THE  SYSTEM 

4 

26  C 

4 

27  C 

A 

-  A  REAL  2-D  ARRAY  OF  SIZE  NDIM*N  HOLDING  THE 

4 

28  C 

ORIGINAL  MATRIX.  THIS  ARRAY  IS  NOT  ALTERED  BY 

4 

29  C 

IMPRV . 

4 

30  C 

4 

31  C 

LU 

-  A  REAL  2-D  ARRAY  OF  SIZE  NDIM*N  HOLDING  THE 

4 

32  C 

LU  DECOMPOSITON  OF  THE  ORIGINAL  MATRIX.  THIS 

4 

33  C 

ARRAY  IS  NOT  ALTERED  BY  IMPRV 

4 

34  C 

4 

35  C 

NPIV 

-  A  INTEGER  VECTOR  OF  SIZE  N  HOLDING  THE  PIVOT 

4 

36  C 

INFORMATION  FROM  THE  ELIMINATION  STEP. 

4 

37  C 

4 

38  C 

B 

-  A  REAL  VECTOR  OF  SIZE  N  HOLDING  RIGHT  HAND 

4 

39  C 

SIDE  OF  THE  ROIGINAL  SYSTEM  TO  BE  SOLVED.  THIS 

4 

40  C 

ARRAY  IS  NOT  ALTERED  BY  IMPRV. 

4 

41  C 

4 

42  C 

X 

-  A  REAL  VECTOR  OF  SIZE  N  HOLDING  THE  INITAL 

4 

43  C 

APPROXIMATE  SOLUTION.  ON  RETURN  THIS  VECTOR 

4 

44  C 

WILL  CONTAIN  THE  IMPROVED  APPROXIMATION 

4 

45  C 

4 

46  C 

Z 

-  A  REAL  VECTOR  OF  SIZE  N  WHICH  IS  DECLARED  BUT 

4 

47  C 

NOT  INITIALIZED  BY  THE  CALLING  PROGRAM.  ON 

4 

48  C 

RETURN  Z  WILL  CONTAIN  THE  CORRECTIONS  TO  THE 

4 

49  C 

GIVEN  APPROXIMATE  SOLUTION 

4 

50  C 

4 

51  C 

R 

-  A  REAL  VECTOR  OF  SIZE  N  WHICH  IS  DECLARED  BUT 

4 

52  C 

NOT  INITIALIZED  BY  THE  CALLING  PROGRAM.  ON 

4 

53  C 

RETURN  "R"  WILL  CONTAIN  THE  RESIDUAL  R=B-AX 

4 

54  C 

WHERE  X  IS  THE  INITIAL  APPROXIMATE  SOLUTION. 

4 

55  C 

DOUBLE  PRECISION  ARITHMETIC  IS  USED  IN 

4 

56  C 

COMPUTING  "R" 

4 

57  C 

4 

58  C 

4 

59 

REAL*8  AA, 

, XX , SUM 

4 

60 

DIMENSION 

A(NDIM,NDIM) ,Z(NDIM) ,R(NDIM) .NPIV(NDIM) 

4 

61 

DIMENSION 

B(NDIM) ,X(NDIM) ,LU(NDIM,NDIMl 

4 

62  C 

4 

63  C 

CALCULATE  THE 

RESIDUALS 

-  32  - 


4 

64  C 

4 

65 

DO  20  1=1 ,N 

4 

66 

SUM  =  B(I) 

4 

67 

jDO  10  J=  1 ,  N 

4 

68 

AA=A( I , J) 

4 

69 

XX=X(J) 

4 

70 

SUM=SUM-AA*XX 

4 

71  10 

iCONTINUE 

4 

72 

:R(I)  =  SUM 

4 

73  20 

CONTINUE 

4 

74  C 

4  75  C  THE  RESIDUAL  SYSTEM  IS  SOLVED 

4  76  C  i 

4  77  CALL  SOLVE (NDIM,N,LU , NPIV , R , Z) 

4  78  C 

4  79  C  THE  IMPROVED  APPROXIMATION  IS  COMPUTED 

4  80  C  j 

4  81  DO  30  1=1, N 

4  82  X(I)=X(I)+Z(I) 

4  83  30  CONTINUE 

4  84  C  i 

4  85  RETURN 

4  86  END 


i 

i 
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